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ABSTRACT 

We study the imprints of boxy/peanut structures on the 2D line-of-sight kinematics 
of simulated disk galaxies. The models under study belong to a family with varying 
initial gas fraction and halo triaxiality, plus few other control runs with different struc¬ 
tural parameters; the kinematic information was extracted using the Voronoi-binning 
technique and parametrised up to the fourth order of a Gauss-Hermite series. Building 
on a previous work for the long-slit case, we investigate the 2D kinematic behaviour 
in the edge-on projection as a function of the boxy/peanut strength and position an¬ 
gle; we find that for the strongest structures the highest moments show characteristic 
features away from the midplane in a range of position angles. We also discuss the 
masking effect of a classical bulge and the ambiguity in discriminating kinematically 
this spherically-symmetric component from a boxy/peanut bulge seen end-on. Regard¬ 
ing the face-on case, we extend existing results to encompass the effect of a second 
buckling and find that this phenomenon spurs an additional set of even deeper min¬ 
ima in the fourth moment. Finally, we show how the results evolve when inclining the 
disk away from perfectly edge-on and face-on. The behaviour of stars born during the 
course of the simulations is discussed and confronted to that of the pre-existing disk. 
The general aim of our study is providing a handle to identify boxy/peanut structures 
and their properties in latest generation IFU observations of nearby disk galaxies. 

Key words: methods: numerical - galaxy: evolution - galaxy: kinematics and dy¬ 
namics. 


1 INTRODUCTION 


Boxy/Peanut bulges (hereafter B/Ps) are characteristic 
structural features visible in a large fraction of inclined disk 
galaxies (for an iconic example, see ESO 597-36). They are 
now believed to arise from vertical instabilities of the xi 
orbital family - elongated, planar orbits sustaining bars 
iContopoulos &: Papayannopoulos|1980 Athanassoula et al. 


1983 Pfenniger|1984| [Patsis, Skokos, fc Athanassoula|2002| 

Athanassoula|2005b etc.) - although 

competing theories on 
3r 60% of disk galaxies 
1101, while the fraction 

their origin existed until recently. Ov 
are found to host a bar (Buta et al.|2C 

of B/Ps is estimated around 45% l 

Liitticke, Dettmar, & 

Pohlen 2000a Yoshino fc Yamauchi 

20141. The two struc- 


tures are rarely seen together, as the inclination exposing 


one tends to conceal the other (but see [Athanassoula fc 
|Beaton| |2006| [Erwin fc Debattista| [2013[ and references 
therein); that their occnrrence in disk galaxies is found to 
be broadly compatible provides in itself a “statistical” piece 
of evidence suggesting a common origin (de Souza & Dos 
Anjos[[l987 Shaw[[l9F7 1. The first numerical simulations of 
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isolated disk galaxies had already witnessed B/P formation 


from vertical thickening of the inner regions of bars (Combes 


fc Sanders|[1981| [Combes et al.[[l990| [Pfenniger fc Friedli 
1991[ Raha et al. 19911 albeit in idealised settings. Later 


on, the density profiles of bar-driven, simulated B/Ps were 
found compatible with observed light profiles along different 


cuts I Athanassoula & Misiriotis 2002 Liitticke, Dettmar, 


fc Polilen[[2000b Aronica et al.[[2003 Athanassoula[[2005b I; 
similarly, unsharp-masked images revealed akin morpho¬ 
logical features (Bureau et al. 20061. On top of these 


photometrical studies, successful comparisons between the 
line-of-sight kinematics of simulated and observed B/Ps 
further disclosed the tie between these structures and bars 
([Kuijken fc Merrifield 1995 Athanassoula fc Bureau 1999 


Bureau fc Freenian|1999||Merrifield fc Kuijken|1999 Chung 


fc Bureau J200^ Kormendy fc Kennicutt] 2004 Bureau & 


Athanassoula[[2005[ hereafter BA05; Debattieta/et^alJ|200& 

hereafter D05; [Mendez- Abreu et al.[[2008| [Mendez-Abreu[ 
et al.|2014t 


The line-of-sight velocity distribution hereafter 

LOSVD) contains all the accessible information in an object. 
Indeed, it directly descends from the system’s distribution 
function / after the insurmountable ignorance on the distri- 
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bution of stars along the line of sight, and on their velocities 
across it, are taken into account: 

Ji’{x±,v\\) = J dx|| J /(x,v)d^vx. (1) 

An approximate description of this function is provided by 
its moments up to some order n. The 0-th, 1st and 2nd-order 
moments correspond to, respectively, the surface brightness 
/i, the mean line-of-sight velocity V and velocity dispersion 
a: 

m(xx) =/^(xx,V||)dv||, (2) 

y(xx) = - [ -^(xx, V||)v||dv|j, (3) 

M J 

cr(xx) = ■^(xu,V||)(v|| -y)^dv||, (4) 

and the higher moments follow suit. 


Observationally the LOSVD is obtained from the 
broadening of spectral lines, given an educated assumption 
on the emitted spectrum and accurate calibration of instru¬ 
mental effects (Kuijken & Merriheld 1993 and references 
therein; Binney &: Merrifield|1998 1. The kinematic informa¬ 
tion it contains is generally quantified in terms of a limited 
set of parameters. These mirror the moments mentioned 
above and control the shape of an ad-hoc analytic function 
which is fitted to the observed curve. 

A Gaussian shape was initially regarded as a good ap¬ 


proximation to observed LOSVDs (Sargent et al. 1977 
|Tonry fc Davis] |1979| ) and this implied a description of 
the kinematics in terms of /r, V and a. When significantly 
non-Gaussian profiles started being observed in the core of 
elliptical galaxies (e.g. |Franx fc Illingworth||1988| ), or even 
in numerical models (e.g. Dejonghe 19871, the Gaussian 
assumption was dropped in favour of a more appropriate 
shape (|Rix fc White] |1992[ |Gerhard| |1993| [Kuijken Sz\ 


Merriheld 1993 Zhao & Prada 19961. The model which is 


now universally adopted was presented in [van der Marelj 
jfc Fran^ ( |1993| ); it describes the LOSVD as a sum of 
orthogonal functions - the Gauss-Hermite series - where 
the Gaussian solution features as a low-order approximation 
to the result. A satisfactory and observationally feasible 
description of the LOSVD in the new parametrisation can 
be obtained by truncating the expansion at the fourth 
order, resulting in two extra parameters besides /r, V and 
(j; these are related to the skewness and kurtosis of the 
LOSVD and are referred to as /13 and /14, respectively. A 
distribution with a more prominent high (low) velocity tail 
will have positive (negative) /13; if its central peak is more 
(less) prominent than for a Gaussian, it will be described 
by positive (negative) /14. 

In summary, we currently reduce all the accessible informa¬ 
tion of a galaxy to the behaviour of fi, V, a, hi,, /14 at as 
many locations as possible on its projected image. 


In this spirit, BA05 studied the kinematic signatures 
of simulated B/P bulges to add yet another piece of infor¬ 
mation that could point to their origin. Their analysis was 
based on collisionless simulations similar to those presented 
in Athanassoula & Misiriotis (2002 1 and focussed on the 


major-axiQ behaviour of highly inclined systems, where the 
B/P shape is fully visible; their results were therefore rel¬ 
evant for comparison to ID, “long-slit” kinematic studies 
such as those of Chung & Bureau (20041. In this paper we 
make use of a newer set of dynamical simulations featuring 
star formation to extend the analysis of BA05 to 2D and dif¬ 
ferent structural components. The purpose is to present the 
results in the form of kinematic maps similar to those pro¬ 
vided by ever more sophisticated and widespread Integral- 
Field observations (see, e.g., Seidel et al.|2014 |. This is par¬ 
ticularly relevant as B/P structures extend above the disk 
and may have considerable 2D cross sections at all disk in¬ 
clinations and viewing angles. In addition, projection effects 
may at times occult their characteristic shape and hinder 


their morphological identihcation (Athanassoula & Beaton 


2006 Erwin & Debattista 20131; in these cases, kinematic 


information is crucial to assess the presence of these struc¬ 
tures in the observed object. Even when the B/P is visible, it 
is difficult to assess its orientation based only on its appear¬ 
ance. Indeed, the relation between its horizontal and vertical 
extent allows for considerable scatter and this prevent ac¬ 
curacies better than 40-50 degrees in the determination of 
the B/P position angle, whence the need for complemen¬ 
tary kinematic information to pinpoint the orientation of 
the structure (Athanassoula & Martinez-Valpuesta, unpub¬ 
lished; [Martinez- Valpuesta fc Athanassoula||2008| and ref¬ 
erences therein). These topics will be addressed in turns, 
according to this scheme: 


Sec.d - we introduce the simulations and the method used 
to obtain the LOSVD; 

Sec. m - edge-on case: we focus on the ID results and com¬ 
pare them to those obtained by BA05; 

Sec.H - edge-on case: we present the 2D maps and dis¬ 
cuss the effect of B/P strength, position angle as well as the 
contribution of bars and B/Ps separately; 

Sec. m - edge-on case: we investigate the masking effect of 
a classical bulge; 

Sec. [H - face-on case: we discuss the signatures of B/Ps 
with different strengths and evolutionary history; 

Sec. IZl - we discuss different system setups; 

Sec. [H - we comment on the effect of disk inclination; 

Sec. m - we summarise the work and lay our conclusions. 


2 METHODS 

In this section we will describe the simulation dataset and 
the technique used to extract the line-of-sight kinematics. 


2.1 Simulations 


We perform the study on simulations of idealised disk 
galaxies, commonly referred to as “dynamical” (see 


Athanassoula 2013 for a review). Specifically, the bulk 


of the runs under investigation belong to the GTR (Gas 


TRiaxial) series presented in Athanassoula, Machado, & 


Rodionov (20131 (hereafter AMR13). This consists of a grid 


of simulations with varying initial gas fraction and halo 


^ By kinematic major axis we mean the cut across the projected 
galaxy which encompasses the highest rotational velocity. 
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triaxiality, developed to test the effect of these two aspects 
on bar formation and growth. We refer to AMR13 for a 
thorough description of the numerical setup and report 
the basic details in Table We study hve of the GTR 
runs; they share all the disk properties (except in that 
the mass contributed by gas vary from 0% to 75% of the 
total), as well as the halo proHle. The latter is a core of 
spherical shape in all but GTR102, where a mild triaxiality 
is present. To this set we add two simulations with a cusped, 
NFW halo ( [Navarro, Frenk, fc White|1996 1, a different disk 
scalelength and initial rotation properties: GGS006 and 
GGS008 (Gas CuSp, see Athanassoula et al.j 20141. They 
act as “control runs” as they lift the homogeneity in the 
structural parameters of the other simulations. Finally, to 
investigate the effect of a classical bulge we add simulation 
CBL005 to the study; here this component is modelled as 
proposed by Dehnen (19931 and we refer to Table for the 
assigned parameter values. 

The number of collisionless disk particle^ varies among 
the runs and is maximum for GCS006, where it reaches 
8 X 10®. The number of gas and stellar particles not only 
changes according to the simulation, but also evolves from 
snapshot to snapshot; GTR116 contains the largest number 
of gas particles at the initial conditions (7.5 x 10®) and the 
largest number of stars at the end (1.38 x 10®). Starting 
from these numbers and following the procedure adopted in 


Athanassoula (2005bI, the particles in both the stellar and 


disk components have been increased by a factor of forty 
before the extraction of the kinematics was performed. This 
was achieved on one side by replicating particles according 
to the four-fold symmetry of the system, a procedure that 
brings a factor four in particle boost. The rest was provided 
by the stacking of ten adjacent snapshots (separated by 
5 Myr in time) after aligning the bar symmetry axes. 
In order to check whether the above procedure altered 
the results, we performed the kinematic analysis also 
on the original snapshots of GTRlOl and GTR116; we 
found no significant difference, except for an increased noise. 


Fig. [^andj^show, respectively, the evolution of bar and 
B/P strength for the GTR runs. 

The former is evaluated by taking the maximum of the rel¬ 
ative amplitude of the m = 2 Fourier components, as in 
AMR13. Here we show the behaviour of the disk component 
only, as this is the only quantity we will explicitly refer to in 
what follows. Indeed, we will primarily discuss disk results 
at relevant times selected according to this and the follow¬ 
ing plot and comment on the stellar behaviour separately. 
The time evolution of the bar strength evaluated for the disk 
and stellar component together can be appreciated in £g. 7 
of AMR13. 


As for the B/P strength, several definitions have been 
proposed and found to reasonably agree with each other 
(Martinez-Valpuesta & Athanassoula 20081. Here we eval¬ 
uate this quantity as the maximum value reached by the 
median of the \z\ coordinate as a function of x for an edge- 


^ Throughout the paper we will use the convention whereby 
“disk” particles are those present, in fixed number, from the start 
until the end of the simulation; “star” particles are, instead, those 
forming progressively in time out of the gaseous component. 


Bar strength vs. time 



Figure 1. Evolution of the bar strength for the disk component 
of the GTR runs under study. This is quantified as the maximum 
of the relative amplitude of the m = 2 Fourier components. 


on view. Note that this measurement is not sensitive to the 
position angle of the B/P because it reflects only its ver¬ 
tical extent. As will be discussed in Sec. we considered 
other percentiles than the median and reached compatible 
results; in general, all the method used compare very well 
with the visual assessment of the maximum vertical extent 
of the structure. In order to allow for comparisons among the 
curves, we normalised this quantity by the vertical extent of 
the outer disk; this was computed as the minimum median 
value of the vertical extent of the particles outside the B/P 
structure and within 15 kpc from the centre. Again, as for 
the bar strength, we show the results for the disk particles 
only. The stellar component does not have a radically dif¬ 
ferent behaviour, however; it starts out with a much lower 
vertical extent than the disk, as expected, and eventually 
catches up with it at late times. 

Based on time evolution of bar and B/P strength, we define 
a strong-B/P and a moderate-B/P group of simulations. The 
former contains simulations GTRlOl, GTR102, GCS006 and 
GTR106 - i.e. those undergoing one or more bucklings, sud¬ 
den drops in bar strengths accompanying the formation of 


strong B/Ps (Raha et al. 

1991 

Pfenniger & Friedli 

1991 

Martinez-Valpuesta & Shlosman|2004||Athanassoula|2005a 

Martinez-Valpuesta et al.|2006 

Debattista et al.|2006 

. The 


latter concerns, instead, the runs presenting a more grad¬ 
ual bar and B/P growth (GBL005, GTR116, GTRlll and 
GCS008). 


With the exception of Sec. we will focus on perfectly 
edge-on and face-on projections. These correspond to a disk 
inclination i = 90 and i = 0, respectively. In both cases we 
rotate the system such that the longer side of the bar-B/P 
runs parallel to the x axis, while the short side to the y 
axis. In the edge-on case we distinguish (i) the side-on view 
{xz plane, position angle 0 degrees), where the long side of 
the bar-B/P is fully visible (ii) the end-on view [yz plane, 
position angle 90 degrees) where, conversely, the short side 
of the bar-B/P becomes exposed. Within this convention, 
by “small position angles” we will therefore refer to close- 
to-side-on views, while by “large position angles” to close- 
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Name 

Halo 

[7 in kpc] 


f^halo 

[ 1 O 1 «M 0 ] 

Disk 

[h, zq in kpc] 

^disk 

[IOIOM 0 ] 

Gas fraction 

Bulge 

[ 7 , a in kpc] 

^bulge 

[IOIOM 0 ] 

GTRlOl 

core, 7 = 

1.5, sphr 

25 

3,0.6 

5 

0 

- 

- 

GTR106 

core, 7 = 

1.5, sphr 

25 

3,0.6 

5 

0.25 

- 

- 

GTRlll 

core, 7 = 

1.5, sphr 

25 

3,0.6 

5 

0.5 

- 

- 

GTR116 

core, 7 = 

1.5, sphr 

25 

3,0.6 

5 

0.75 

- 

- 

GTR102 

core, 7 = 

1.5, triax 

25 

3,0.6 

5 

0 

- 

- 

GCS006 

nfw, C 200 

= 8.3, sphr 

19 

4,0.6 

5 

0 

- 

- 

GCS008 

nfw, C 200 

= 7.5, sphr 

24 

4,0.6 

5 

0.4 

- 

- 

GBL005 

core, 7 = 

15, sphr 

44 

2.5,0.4 

6.4 

0 

1.5,0.9 

1.8 


Table 1. Structural features of the simulations under study. For a cored halo 7 represents the core radius, while in the NFW case C200 
represents the concentration. For the latter we report the total mass instead of the typical M 20 O 1 a-s the distribution is truncated at 
around half of the corresponding virial radius. The tags “sphr” and “triax” (for spherical and triaxial) describe the halo shape. The 
behaviour of the disk component is described by an exponential function in radius and by a sech^ in the vertical direction; h and zq 
represent the scalelength and scaleheight, respectively. The gas fraction expresses how much of the initial disk mass is in the form of 
potentially star-forming, SPH particles. The classical bulge is modelled as described in |Dehnen| ( [1993^ ; 7 controls the slope of the inner 
profile and a is its scalelength. 


Peanut strength vs. time 



Figure 2. Evolution of the B/P strength for the disk component 
of the GTR runs under study. This is quantified as the maximum 
of the median \z\ value computed along the x direction when the 
B/P structure is viewed side-on. To allow comparisons between 
the runs, this measurement is normalised to the median \z\ value 
of the outer disk. 


to-end-on views. Finally, all models rotate counterclockwise 
when seen face-on. 


2.2 Extracting kinematics 


Once the simulated galaxy is placed at the desired orienta¬ 
tion, the first step towards deriving its kinematic behaviour 
consists in defining the LOSVD. This entails a decision on 
how to group particles together to give rise to a single local 
velocity distribution, a procedure known as “binning” and 
which is crucial for deriving reliable kinematic information 
from the object. Here we will use the now widespread 


method developed by Cappellari & Copin (20031 and 


known as “Voronoi-binning technique”. The starting points 
are a given held-of-view (hereafter FOV) where the signal 
is pixelated and a requirement on the minimum desired 
signal-to-noise (hereafter S/N). When the method is applied 
to simulations, the signal corresponds to the number of 


particles in each pixel; similarly, a requirement on S/N 
translates into a requirement on the number of particles 
under the assumption of Poissonian noise {S/N = FS)- 
The desired S/N level decides how many particles one 
needs to group in order to obtain a LOSVD which is 
robust enough to allow for the desired accuracy in the 
derived kinematic quantities. However, the downside of 
requiring a high accuracy is the loss in spatial resolution; 
consequently, values in excess of S/N = 50 are hardly 
ever used as this is commonly considered sufficient to 
obtain reliable information on the higher moments /13 and 
/14 ( [van der Marel fc Fran:j^|1993[ ). Given these elements, 
the Voronoi method optimally tessellates the FOV, i.e. it 
preserves the maximum spatial resolution whilst fulfilling 
the S/N criterion for each of the resulting cells. Each of 
these regions will give rise to a LOSVD which is then 
htted to the Gauss-Hermite series introduced in Sec. 
as a consequence, each of the Voronoi bins will have an 
associated V, a, ha, ^4 - which is what is finally shown in 
the 2D maps. The units used here are km/s for the first and 
second moment, while ha and h ,4 are dimensionless. As for 
the surface brightness, or “SB”, we will just express it as 
the logarithm of the total mass (in Mq) contained in each 
pixel. 


The procedure we used to perform the gridding of the 
data down to the extraction of the kinematic maps stems 
from that developed by Brown et al. (20131 and differs from 
the original only in its optimisation. We adapted it to run on 
ninety snapshots for each component and eight B/P view¬ 
ing angles of each simulation (resulting in between 720 and 
1440 sets of kinematic maps per simulation); because in 
some cases the number of particles undergoing the analy¬ 
sis exceeded fifty million, both wallclock time and memory 
requirements had to be downsized accordingly. 

We kept the FOV and pixel size uniform among the strong 
and moderate-B/P groups: in the former they are set to 
20 x 10 kpc^ and 0 . 1 x 0.1 kpc^, respectively, while in the latter 
to 16x8 kpc^ and 0.07x0.07 kpc^. The FOVs are purposely 
large to encompass the entire B/P region in the side-on view 
as well as part of the outer disk; indeed, the aim here is not 
to reproduce a specific observation but rather to identify 
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B/P imprints thoroughly. The pixel sidelength is 40 to 100% 
larger than the adopted softening length e = 0.05 kpc. As an 
example, our choices compare well to the FOV and spatial 
resolution of the MUShlflinstrument in Wide-Field-Mode for 
an observation of a nearby object. Indeed, at a distance of 
48 and 38 Mpc the median resolution of 0.46 arcsec (valid 
at 0.75pm and with adaptive optics) translates into sizes of 
0.107 and 0.084 kpc, respectively - in broad agreement with 
the pixel values used here. 

Finally, for both sets of simulations we require S/N = 50 
although we also re-ran cases with different S/N levels as 
a sanity check. As mentioned above, the higher the S/N 
the more accurate is the evaluation of the kinematic mo¬ 
ments - especially the high-order ones, which are sensitive 
to the behaviour of the tails of the distribution. This holds, 
of course, under the assumption that the Gauss-Hermite se¬ 
ries is a good description of galactic LOSVDs. One should 
also consider that as they increase in size the Voronoi bins 
will start encompassing regions that may have intrinsically 
different kinematic features - resulting in a blending of the 
signatures. 

In Fig. we show, as an example, the variation of the /14 
maps for simulation GTRlOl at t = 8 Gyr when different 
S/N values are adopted. We only show /14 as it is the quan¬ 
tity whose behaviour is most affected by changes in S/N. 
Together with it, we show the errors provided by the fit¬ 
ting procedure. A blending effect is evident, as different re¬ 
gions progressively merge into the same bin. However, the 
morphology of the maps is qualitatively unchanged as the 
regions of minima/maxima maintain their identity and rel¬ 
ative positions. The errors for the case S/JV = 50 can be 
comparatively large in a handful of bins, but never enough 
to overturn the interpretation of the global /i 4 behaviour. As 
expected, they reduce in importance for higher S/N values. 
The bins presenting a formal null error, and correspondingly 
the most negative /14 values, present markedly non-Gaussian 
LOSVDs; as a consequence, the fitting procedure hits the 
lower bound set for /14 and no formal error is returned. This 
problem is found to affect mainly h 4 , while ha hardly ever 
attempts taking on extreme values. Although the h 4 results 
eventually assigned to the affected bins should not be taken 
at face value, the qualitative behaviour they suggest finds 
support in that of the neighbouring, well-behaved cells. 
Overall, we find that setting S/N = 50 provides converged 
kinematic information - the more so the lower the order 
of the moment under consideration. We warn against tak¬ 
ing extreme /14 values literally, as we often find these being 
related to non-Gaussian LOSVDs, and rather focus on the 
qualitative behaviour they point to. 


3 EDGE-ON VIEW: ID RESULTS 

In this and the next two sections we investigate the kine¬ 
matic properties in the edge-on view. The reference work 
for this is BA05. 

^ The Multi Unit Spectroscopic Explorer (MUSE) is 
an integral-field spectrograph operating at visible wave¬ 
lengths, mounted at the beginning of 2014 on the 
UT4 telescope of the ESO VLT facility in Paranal: 
http: / / www.eso.org/sci / facilities/develop/instruments / muse.html 


In Fig. we show the major-axis kinematics for simulation 
GTRlOl at t = 8 Gyr. The different lines show the varia¬ 
tions with the B/P position angle. The main B/P-induced 
features discussed by BA05 can also be appreciated in this 
plot. They are: 


• a quasi-exponential light profile, with a steep descent 
and a plateau at large radii; 

• a double-hump velocity profile characterised by the 
presence of a local minimum/maximum; 

• a broad, flat peak in velocity dispersion, accompanied 
by secondary maxima farther out; 

• a correlation between ha and V extending till the end 
of the B/P structure; 

• a fiat central minimum in /14 followed by a rise and 
decline at larger projected radii. 


These global kinematic properties of B/P structures ul¬ 
timately arise from the behaviour of the orbital families they 
originate from, most notably the Xi group. In particular, 
the differences among the curves in each panel can be un¬ 
derstood in terms of the viewing-angle dependence of the 
various classes of orbits ( Bureau fc Athanassoula]|1999 1. In 
the end-on projection, where the line-of-sight kinematics is 
dominated by the Xi family, the aforementioned features 
reach their maximum strength. 

BA05 remark that not only are these kinematic features 
more prominent for larger position angles, but also for 
stronger B/P structures. Fig. compares the end-on re¬ 
sults from three simulations with decreasing B/P strength 
at t = 8 Gyr: GTRlOl, GTR106 and GTR116. For the 
last two we show the behaviour of the stellar component 
alongside that of the disk. Star formation induces quite dra¬ 
matic changes in the density profiles, which become more 
pointed; the relative importance of the contribution by the 
two components depends on the initial gas fraction in the 
simulation. The double-hump structure in velocity is weaker 
for GTR106 and GTR116, and in the latter the local mini¬ 
mum/maximum is hardly visible. The behaviour in veloc¬ 
ity dispersion reflects the different projected sizes of the 
B/P structures, but the plateau and secondary maxima are 
present in all the three cases. Similar considerations hold for 
the extent of the hs-V correlation region, although we note 
that in GTR116 hs is also slightly reduced in magnitude. As 
for /i 4 , the central minimum and the subsequent rise become 
less important for weaker B/Ps; same holds for the minima 
farther out, deeper and wider for stronger B/Ps and whose 
position reflects the projected size of the structure. 

Aside from the different contributions to the projected den¬ 
sity, the disk and stellar component have similar behaviours 
along the kinematic major-axis. The only notable difference 
is found in the a and /14 profiles for GTR106: the stellar 
component presents a central dip in velocity dispersion, re¬ 
flected in an ^4 excess. A similar feature was found by BA05 
in their strong-B/P, collisionless simulation at large viewing 
angle. This was explained as a consequence of the xi orbit 
becoming more circular towards the centre, and this may 
indeed explain the slight presence of this feature in GTRlOl 
at large position angles (Fig. [^. However, in GTR106 the 
effect concerns essentially only the stellar component and a 
dissipational process must be invoked. As it turns out, the 
sigma dip is induced by young stars (born at t > 2 Gyr) and 
reflects the temperature profile for the gas at late times; a 
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Figure 3. End-on view of simulation GTRlOl at t = 8 Gyr. The top rows shows the variation of the maps for the disk component 
when adopting S/N = 50,100,150,200 (from left to right). For each of the cases, the bottom rows shows the errors provided by the 
fitting procedure in each Voronoi bin. 



Figure 4. Major-axis kinematics of the disk component from 
simulation GTRlOl in an edge-on projection. The simulation is 
taken at t = 8 Gyr, i.e. well into the secular stage of B/P evolu¬ 
tion. The panels show the mean behaviour of the five kinematic 
quantities along a 700-pc-wide slit centred on the z = 0 plane. 
Each line corresponds to a different position angle, from 0 (side- 
on view) to 90 (end-on view) in steps of fifteen degrees. The grey 
line in the top panel marks the disk’s original profile for reference. 


Figure 5. Major-axis kinematics for simulations GTRlOl, 
GTR106 and GTR116 in an end-on view. The simulations are 
ordered by decreasing B/P strength at t = 8 Gyr. For GTR106 
and GTR116 we show separately the disk and stellar components 
(solid and dotted lines, respectively). 


similar, albeit less striking behaviour characterises GTRlll 
and GCS008 (not shown) while it is absent in GTR116. Dif¬ 
ferences in the behaviour of the disk and stellar components 
are more evident in 2D; we will discuss this aspect in Sec. |4.4| 

Our ID results agree with the picture outlined by BA05. 
The strength of the B/P structure and its position angle de¬ 
termine how striking the projected kinematic properties will 
be. The degeneracy between these two quantities is limited, 
though. From Fig. and it is not obvious to distinguish 
GTR106 end-on from GTRlOl at some smaller position an¬ 
gles and indeed some ambiguity in disentangling the effect 


of the two quantities is present. However, the behaviour of 
GTR116 is such to prevent the identification with a much 
stronger B/P structure seen at a smaller position angle. The 
magnitude reached by /13 in particular, points in itself to 
(close to) an end-on view. 

In summary, B/P strength and position angle influence the 
line-of-sight kinematics in ways that are not entirely over¬ 
lapping. The highest moments of the LOSVD, in particular, 
evolve differently with increasing B/P strength and position 
angle; this will be clearer from the results of the following 
section. 
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4 EDGE-ON VIEW: 2D KINEMATIC MAPS 

We will now show how B/P kinematic signatures evolve 
away from the kinematic major axis. As in the previous sec¬ 
tion, we will start by investigating the variations induced by 
the B/P position angle for individual simulations; we will 
then contrast all the simulations at a given position angle to 
assess the effect of B/P strength. This comparison is done at 
at t = 6 Gyr, but the choice of time is not crucial to assess 
the relative differences between the simulations once these 
have entered the slow phase of B/P secular evolution. The 
reason for this particular choice will be clarified in Sec. 
where we will discuss the masking effect of a classical bulge. 
Until then we will not explicitly discuss simulation CBL005. 
We will conclude the section with a comparison between bar 
and B/P-induced kinematic signatures, as well as a discus¬ 
sion on the behaviour of stars. 

As a reference, we show in Fig. the appearance of the 
maps a.t t — 0.05 Gyr - i.e. at a moment where all the runs 
are still axisymmetric. Comparison against this panel should 
ease the isolation of the kinematic features which are gen¬ 
uinely induced by bar/peanut structures. We will comment 
on some of these aspects at the end of the section. 


4.1 The effect of position angle 

Fig.[7] and show, respectively, simulations GTRlOl and 
GTR116 at different position angles. We decided to show 
these two runs as they are quite representative of the strong 
and moderate-B/P groups. 

As far as the velocity behaviour is concerned, we recover 
some of the features pointed out by [Athanassoula fc] 


Misiriotis (2002 1 for simulations of similar scope as those 


used here. As expected we find some degree of cylindrical 
rotation in the B/P region, suggesting these structures 
move coherently to some extent. As already noticed by 


Combes et al. ( 

1990|), Athanassoula & Misiriotis (20021 

and reported by 

Williams et al.| |[2011[), when going from a 


side-on to an end-on view the spacing between a given set of 
isovelocity lines progressively shrinks towards the centre. In 
some cases this is accompanied by a decrease in the degree 
of cylindrical rotation - i.e. the isovelocity lines evolve from 
a parallel to a splayed set. This is the case of GTRlOl, for 
example, while for GTR116 the effect is much weaker. The 
idea of cylindrical rotation as an ubiquitous feature of B/P 
bulges has been challenged from multiple angles: |Athana^ 
|soula fc Misiriotis] ( |2002| showed how models with different 
disk and halo contributions to the total velocity curve may 
result into B/P structures with different rotation properties 
at different position angles; [Williams et al.| ( |2011| ) analyse a 
sample of five edge-on disks with B/P structures showing a 
considerable range of behaviours as far as projected velocity 
is concerned. In our analysis we register a tendency for 
cylindrical rotation to weaken when going from side-on to 
end-on views, although the importance of this effect changes 
from case to case without a clear trend with B/P strength. 
Similarly, at a given position angle the degree of cylindrical 
rotation of the different runs does not follow the sequence 
in B/P strengtlQ Overall, even though we attempted a 


The reasons for this lack of clear trend can be several: the 
quantification of cylindrical rotation comes with some arbitrari- 


quantification of the evolution of the isovelocity lines as 
a possible mean to identify the position angle of a B/P 
structure, we find that this diagnostic does not significantly 
improve the results discussed below and will not indulge on 
this aspect further. 

Moving on to higher moments, a common feature of the 
GTRlOl and GTR116 sequences is the gradual evolution of 
the a maps, reflecting the larger amount of orbital families 
and cuspier orbital shapes contributing to the LOSVD 


at large position angles (Athanassoula 1992 Bureau & 
Athanassoula|19OT I. This effect was also captured by the ID 


analysis and indeed the major-axis behaviour summarises 
adequately the global evolution of the second moment. 
When considering /is and /i 4 , however, the evolution is more 
dramatic for stronger B/Ps and involves interesting features 
offset from the major axis. In this case, the kinematic 
maps change considerably from the side-on to the end-on 
projection in both the amplitude of the moments and their 
morphology; starting from a bland appearance in the side-on 
view, at intermediate position angles the higher moments 
develop elongated “wings” in spatial correspondence with 
the projected B/P edges. These tend to disappear from the 
/is maps when approaching the end-on view and instead 
evolve into a strong X shape in /14. For moderate B/Ps 
the changes are considerably reduced and are essentially in 
the magnitude of the moments more than in the general 
morphology. The region of h^-V correlation increase in 
size and magnitude, maintain the overall appearance of 
homogeneous blobs extending just below and above the 
disk plane. The /14 minima grow deeper and become better 
defined spatially, but no outstanding signatures arise at any 
position angles. 


To further stress how the position-angle dependence of 
the higher moments varies with the importance of the B/P 
structure, Fig. [^summarises the behaviour of the / 13 -U cor¬ 
relation for all the simulations under consideration. The lines 
represent the number of pixels at a given abscissa value with 
I/ 13 I > 0.15, i.e they follow the evolution of the regions with 
the largest /i 3 amplitudes. The size of the correlation blobs 
increases notably from the side-on to the end-on view; con¬ 
sidering that all the simulations present similarly feature¬ 
less side-on maps (as will be shown below), the evolution is 
considerably more important at larger B/P-strength values. 
The differences between the runs diminish as we lower the 
threshold on |/i 3 |, but the considerations on the size of the 
affected regions and the dependence on position angles and 
B/P strength still qualitatively hold. Similar results are ob¬ 
tained for /14 even though the shape of the curves becomes 
more irregular due to the morphology of the minima in the 
maps; we do not show these results for conciseness. 


ness (e.g. the choice of the area over which it is evaluated) and 
the B/P strength itself is a more complex quantity than can be 
described by one single number; moreover, the projected velocity 
will carry information on the orbital structure in each bar and on 
various properties of the periodic orbits of the main families. 
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Figure 6. Edge-on view of the disk component at t = 0.05 Gyr, before non-axisymmetric structures develop in the simulations. Note 
that the vertical size of each map varies from simulation to simulation and has exceptionally been reduced with respect to the values 
introduced in Sec. |2.2| to avoid displaying empty regions. The magenta frame emphasises the main GTR family. The units adopted in 
this and subsequent panels are as in Fig. (see also Sec. |2.2| |. 


Effect of position angle on the h 3 -V correlation - t = 6 Gyr 



X -> y [kpc] 


Figure 9. Number of pixels with j/isj > 0.15 as a function of the 
abscissa value (e.g. x in the side-on case) for a number of position 
angles and for each of the runs. The curves are normalised to the 
total number of pixels along the z direction. Each panel represents 
a simulation and each of the coloured curves a different position 
angle, as stated in the legend. In this respect, the GTRlOl and 
GTR116 panels are derived from the 2D results shown in Fig. 
and[8] 


4.2 The effect of B/P strength 

Fig. and m show, respectively, the side-on and end-on 
views of all the simulations under consideration, ordered by 
decreasing B/P strength. The appearance of the kinematic 
maps for the side-on case does not change dramatically 
among the runs (except for CBL005, a case which will be 
discussed elsewhere). The exception is the a map, whose 
“richness” (i.e. peak value and extent of intermediate-cr 
area) is a monotonic function of B/P strength. This is 
better appreciated when focussing on the main GTR family, 
whose panels are framed in magenta; indeed, GCS006 
slightly breaks the sequence due to a slower rotation at 
the initial conditions. An hs-V correlation is present in all 
cases in the inner parts, but the amplitude reached by /is 
is rather limited. The exception is GGS006, although this 
is a transient effect (indeed the isodensity contours suggest 
this simulations may be undergoing some reassessment 
at this time). The behaviour of hi is such that minima 
are present at the end of the B/P shape, while positive 
regions of different shapes dominate the centre. Overall, the 
behaviour along the kinematic major axis discussed in the 
Sec. is quite representative of the 2D picture. 

In the end-on view, the effect of B/P strength is consid¬ 
erable and results in a rather neat sequence of kinematic 
behaviours. This should not surprise given the result of 
Sec. |4.1| whereby the evolution from side-on to end-on 
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Figure 7. Edge-on view of the disk component from simulation GTRlOl at t = 6 Gyr and for different position angles. 


induces more dramatic changes for stronger B/Ps. This 
aspect is exemplified in Fig. which, similarly to Fig. 
summarises the ha behaviour for all simulations at different 
position angles. All runs initially present rather bland 
features and a well-defined sequence starts at position 
angles « 60; the h^-V correlation blobs induced by stronger 
B/P structures start covering larger areas and extend 
farther out from the centre. 

The relative behaviour of the runs in cr and ha for both 
the side-on and end-on views is summarised in Fig. 1 and 
2 from the online material; these show the peak value and 
spatial extent of the main morphological features in the 
maps. 

In summary, we find that 2D maps present a few inter¬ 
esting B/P-related features on top of those already identi¬ 
fied in major-axis studies. These come in the form of wings 
of ha-P correlations as well as of deep, X-shaped hi min¬ 
ima. They are more important for stronger B/Ps and at 
intermediate/large position angles. Even at moderate B/P 
strengths, due to the bin-to-bin variation of the kinematic 
information and mild asymmetries in the B/P features, we 
find that the 2D information allows a better assessment of 


the presence and importance of a B/P. 

Overall we find that for the group of simulations with a 
strong B/P structure the differences generated by the posi¬ 
tion angle are to some extent separable from those induced 
by B/P strength. Having inspected edge-on maps varying 
the position angle in steps of five degrees, we find that a 
rough identification of the structures as being at “small” 
(0 — 30 degrees; bland h^-V correlation and /14 maps), “in¬ 
termediate” (30 — 60 degrees; off-plane features in both /is 
and /14) and “large” (60 — 90 degrees; large /is and /14 am¬ 
plitudes, off-plane /14 features) viewing angles seems feasible 
in principle. 

For the group of simulations with a moderate B/P this is 
true to a lesser extent. We find that the evolution with posi¬ 
tion angle is stronger than the evolution with B/P strength, 
but it is still quite moderate and does not allow to bracket 
the position angle as in the cases discussed above. To compli¬ 
cate matters further, we find that a bland h^-V correlation 
and shallow /14 minima are not exclusively related to the 
presence of bar or B/P structures (see Fig. [^; a region of 
/13-P correlation, in particular, can arise in the central re¬ 
gions of perfectly axisymmetric systems depending on the 
shape of the rotation curve as well as the amount of veloc- 
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Figure 8 . Edge-on view of the disk component from simulation GTR116 at i = 6 Gyr and for different position angles. 


ity scatter allowed around il[^ These pre-existing features 
of the disk may influence the ha behaviour even after the 
formation of weak B/Ps and particularly so at low position 
angles, where the contribution of the B/P orbits to the high- 
velocity tail is smaller. This source of uncertainty, on top of 
the already limited dependence of the kinematic behaviour 
on the viewing angle, makes it is more difficult to estimate 
the position angle of the B/P structure in this case. 

To better appreciate to what extent strength and position 
angle are degenerate in dehning the kinematic behaviour of 
a B/P bulge. Figs. 3 to 7 of the online material show each of 
the five moments separately as a function of the two quan¬ 
tities for all the simulations under consideration. 

Finally, we stress also that only a single model with classical 
bulge is discussed here. While it is unlikely that such axisym- 
metric models can give rise to similar features as those of 
B/P bulges (in particular those clearly associated with the 
B/P edges), we cannot state with certainty that the prop¬ 
erties of CBL005 are representative of all classical bulges. 
This would require generating an extended grid of models, 
something which is beyond the scope of this paper. 

® See also the discussion in Sec. 5 and 7.3 of BA05 


4.3 Individual contribution of bars and B/Ps 

In the secular evolution phase, both the bar and the B/P 
bulge grow stronger in time (see Fig. and Fig. [^. The 
maps evolve accordingly and we wonder if they do so due 
to the growth of one of these two structures preferentially. 
In order to tackle this, we need to know what features 
each of these two structures induce in the kinematics. In 
this section, we inspect the kinematic behaviour of the 
simulations at times of equal bar strength and with or 
without the presence of a B/P. This comparison is only 
feasible for the strong-B/P group of simulations, where - 
due to the buckling phenomenon - the bar strength is not a 
monotonic function of time. We can then identify pre- and 
post-B/P snapshots characterised by identical bar strengths 
and compare their kinematic behaviour as a function of 
position angle. Fig. shows the results for simulation 
GTR102. We picked this run because the behaviour of the 
bar strength is such that two suitable snapshots can be 
found at a large enough distance from the buckling “peak” 
- t — 2 and t — 4.5 Gyr. However, all the strong-B/P 
runs provide very similar pictures as far as this analysis is 
concerned. The main features are listed below. 
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Figure 10. Side-on view of the disk component at t = 6 Gyr. The simulations are ordered by decreasing B/P strength and the magenta 
frame emphasises the main GTR family. 


(7 - Perhaps not surprisingly, at all position angles a in¬ 
creases considerably when a B/P is present. This is due to 
the larger number of orbital shapes populating the struc¬ 
ture. 

hi - The presence of a strong hz-V correlation seems to be 
related to the B/P, more than to the bar structure; in the 
latter case the bland correlation present in the side-on view 
progressively shrinks together with the projected size of the 
bar, to the point that almost no correlation at all is visible 
in the end-on view. The situation is radically different when 
a B/P is present. 

hi - A bar alone induces rather relevant minima in hi es¬ 
pecially at large position angles. However, these are system¬ 
atically stronger when a B/P is present; in addition, a B/P 


induces the four wings of hi minima away from the kine¬ 
matic major axis which are absent in the bar-only case. 

In summary: B/P bulges induce /13 and hi features off the 
kinematic-major axis which, as expected, the bar alone can¬ 
not produce. We find, in addition, that the h^-V correlation 
and hi minima around z = 0 are boosted when a B/P is 
present at a given bar strength. 

4.4 The behaviour of stars 

The stellar counterpart of Fig.[^and[^for the simulations 
with active star formation (GTR106, GTR116, GTRlll, 
GCS008) is given online (Fig. 8 ). Here we will just discuss 
an example from GTR116, the simulation with the most co- 
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Figure 11. End-on view of the disk component at i = 6 Gyr. The simulations are ordered by decreasing B/P strength and the magenta 
frame emphasises the main GTR family. 


pious stellar component. Fig. shows its end-on view at 
t = 10 Gyr. The full disk and stellar populations are pre¬ 
sented in the top two rows, while in the bottom two stars 
belong to different age brackets. Compared to, for example, 
the t = 6 Gyr snapshots (Fig. 8 from the online material), 
the disk and stellar components are not too dissimilar at this 
stage. Note how, for example, the values of a for stars have 
caught up with those for the disk and are only systemati¬ 
cally lower than the latter along the 2 = 0 plane. However, 
stars still carry a somewhat more extreme behaviour in the 
highest moments - especially their amplitude. The even mo¬ 
ments a and /14 are known to be strongly related and indeed 
the presence of positive/negative hi coincides spatially with 
regions of low/high a. This and the lower cr may explain the 


occurrence of prominent regions of positive hi in the stellar 
component. 

In addition, we find that stars born at different times dur¬ 
ing the simulations may present different signatures. The 
behaviour of stars born at t < 0.15 gigayears (i.e. before 
bar/peanut formation) of the start of the simulation is con¬ 
trasted to that of the stars born at t > 2 Gyr (during and af¬ 
ter bar/peanut formation) in the bottom two rows of Fig. |14[ 
each group contains around 15% of the total number of star 
particles at the end of the simulatioij^ The SB maps shows 

® The gas-consumption rate for the GTR runs is visible in fig. 2 
of AMR13; it varies according to the initial gas fraction and for 
GTR116 is particularly high during the early stages of evolution. 
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Figure 13. Edge-on view of the disk component from simulation GTR102 at t = 2 Gyr (no B/P; odd rows) and t = 4.5 Gyr (with B/P; 
even rows), for a selection of position angles. 


already that the two populations have a different spatial dis¬ 
tribution; the youngest stars are preferably found in a very 
cold disk, as suggested also by the behaviour in a. Their cen¬ 
tral line profile is asymmetric, as described by the /13 values 
and considerably more peaked than a Gaussian; indeed, the 
fit attempts to assign /14 values in excess of the fiducial limit. 


5 EDGE-ON VIEW: THE EFFECT OF A 
CLASSICAL-BULGE COMPONENT 

When seen end-on, a B/P bulge loses its characteristic shape 
and acquires a classical-bulge-like morphology. In principle 


the two structures could be distinguished by their light pro¬ 
files, because classical bulges are much more concentrated 
than bars and B/P structures. In practice, however, some 
ambiguity often remains in their identification (see, e.g. 
Gadotti||2009 1 especially as both structures may be present 


at the same time. In this respect, the kinematic behaviour 
can provide a useful piece of information to identify what 
is responsible for the observed morphology. 

In this section we focus on simulation CBL005, a run with 
a moderate classical bulge {Mt^ige/Mdisk+buige ~ 0 . 2 ). 
Because we are interested in its camouflaging effect, unless 
explicitly stated otherwise all the kinematic results shown 
here are for both the disk and bulge component taken 
together. 


At t = 2 Gyr the percentage of disk mass in the form of gas 
has reduced from the initial 75% down to 16%; in the remaining 
8 Gyr this decreases further to 6%, at a much lower pace. 


The B/P structure developed by CBL005 evolves little 
with time and therefore it suffices to show one snapshot from 
this simulation and compare it to the end-on view of all the 
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Figure 14. End-on view of simulation GTR116 at t = 10 Gyr. The first row shows the disk component, while the bottom three different 
stellar samples: all stars, stars born at t < 0.15 Gyr and stars born at t > 2 Gyr, respectively. The side-on counterpart of this panel is 
shown online (Fig. 9). 


Effect of B/P strength on the h 3 -V correlation - t= 6 Gyr 



Figure 12. Number of pixels with jhal > 0.15 as a function of 
the abscissa value (e.g. x in the side-on case) for all the simula¬ 
tions at a number of position angles. The curves are normalised 
to the total number of pixels along the 2 : direction. Each panel 
represents a position angle and each of the coloured curves one of 
the simulations, as stated in the legend. In this respect, the first 
and last panel are derived from the 2D results shown in Fig. |10| 
and llll 

other runs. We picked t = 6 Gyr because it is distant enough 
from the buckling phase of the strong-B/P group of simula¬ 
tions and, at the same time, it is such that the B/P strength 
of CBL005 is somewhat intermediate among all runs. We 
want to see if the presence of a classical bulge can be recog¬ 
nised from the line-of-sight kinematics or, equivalently, to 
what extent the B/P kinematic features are altered due to 
this extra component. 

From Fig. [E where the simulations are ordered by decreas¬ 


ing B/P strength, CBL005 stands out in terms of kinematic 
behaviour. Strong signs of non-cylindrical rotation and large 
velocity dispersion away from z = 0 are not unexpected side 
effects of the presence of a classical component. In terms 
of higher moments, the h^-V correlation is much less pro¬ 
nounced than in most of the runs, while /i 4 takes predom¬ 
inantly positive values in the central regions. If we were to 
order the simulations in terms of decreasing prominence of 
/13 and /14 features, we would probably place CBL005 just 
below GCS008 (see also Fig. |^. Indeed, the two simula¬ 
tions show rather similar high-moment maps; the surface- 
brightness, velocity and cr behaviours are radically different, 
though. In particular, the B/P structure in GCS008 hardly 
sticks out of the disk. 

In Fig. we show the major-axis behaviour for a rep¬ 
resentative subset of the simulations. Alongside GTRlOl, 
GTRlll and GGS008, we show the results of CBL005 with 
and without bulge particles. The shape of the density curves 
is quasi-exponential in all cases except for CBL005 - where a 
central excess is registered. As mentioned at the beginning of 
this section, this is expected. As for the velocity behaviour, 
it is interesting to see how the B/P-driven, double-hump fea¬ 
ture is more evident in CBL005-disk than in CBL005 (and 
perhaps even GTRlll and GCS008); this hints to a mask¬ 
ing effect of the bulge - i.e. when its contribution to the 
LOSVD is added, the B/P-related signatures are damped. 
The simulations with star formation present a more pointed 
a shape, not dissimilar from that of CBL005; GBL005-disk, 
in particular, shows a prominent a peak, and a correspond¬ 
ing hi trough. The reason for this is easily understood when 
looking at the line profiles of the central region. For GBL005 
this is close to Gaussian, (see the small values taken by /is 
and hi) where the central, low-velocity peak is contributed 
by the bulge component and the tails by the disk; when 
removing the bulge the distribution flattens at the centre 
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Effect of a classical bulge; major-axis kinematics 
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Figure 15. Major-axis kinematics for simulations GTRlOl, 
GTRlll, GGS008 and GBL005 in an end-on projection. The 
simulations are taken at t = 6 Gyr, when the B/P strength of 
CBL005 is intermediate between GTRlOl and GTR111/GGS008. 
The panels present the mean behaviour of the five kinematic 
quantities along a 700-pc-wide slit centred on the 2 = 0 plane. The 
results are derived from the disk component except for CBL005, 
where the solid cyan line shows the behaviour of both disk and 
bulge taken together. 


Effect of a classical bulge: offset slit 
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Figure 16. As in Fig. |15| but for a slit centred on 2 = 1.5. 


while maintaining the tails unchanged, leading to a larger 
characteristic dispersion and negative hi. 

From the behaviour of the higher moments one would de¬ 
duce that the B/P structure in CBL005 is weaker than in 
any of the other cases shown, although in reality this is true 
only with respect to GTRlOl. A similar conclusion holds for 
CBL005-disk, albeit to a smaller extent (see the differences 
in V and hs), meaning that the more moderate B/P features 
cannot simply be a consequence of a masking effect enacted 
by the bulge. This is however present and particularly so 
when looking away from the major axis, as shown in Fig. |16| 
Here we inspect the kinematic behaviour of the same simu¬ 
lations along a slit centred on 2 = 1.5. The changes induced 
on the higher moments - especially /13 - by adding/removing 
the bulge are rather striking. In summary, even though the 


presence of a B/P structure in CBL005 can still be recog¬ 
nised from the kinematics, this points to a weaker B/P than 
what really present in the simulation. 

When looking at the morphology of the object, how¬ 
ever, it is evident that a B/P structure cannot be the only 
source of the observed kinematics. As shown in Fig. |17| if a 
B/P were the only responsible for the central protuberance 
then the kinematic maps would look rather different. We 
have selected three snapshots with similar B/P strength 
to CBL005 at t = 6 Gyr from the simulation where this 
was possible. The differences, especially in the highest 
moments, speak by themselves. Note in particular the 
magnitude reached by the h^-V correlation blobs, as well 
as the depth and morphology of hi minima. Fig. |18| and |19| 
show, respectively, the z — 0 and 2 = 1.5 kinematics for 
these runs. The behaviour of hs and hi is rather subdued 
compared to what the CBL005 B/P strength would suggest. 
The velocity dispersion matches well, or even exceeds, those 
of the other runs. The shapes of the density curves are 
similar, except the usual excess for CBL005 at 2 = 0. 

In summary, the kinematic features of a B/P embedded in a 
classical-bulge component are radically different from those 
of a single B/P of the same apparent size. Contributing 
to this effect is (i) an intrinsic modihcation of the disk 
kinematics due to the presence of the bulge (ii) a masking 
effect on the LOSVD. For the case under consideration, the 
presence of a classical bulge component can be guessed by 
(besides the excess of light at the centre): (i) non-cylindrical 
rotation, (ii) weak hs and hi maps compared to what 
expected given the morphology of the structure and the 
high velocity dispersion. 


The question of course is how general this conclusion 
is. The bulge we are considering here is of modest size when 
compared to the typical values found for SDSS galaxies 
{BjT ~ 0.4, Gadotti 20091. It is, however, more promi¬ 
nent than the mean value obtain ed at higher resolut i on for 

(120141 


Laurikainen et al. 


thirty S^G/NIRSOS galaxies by 
(BjT ~ 0.1), who performed the 2D decomposition by in¬ 
cluding a more appropriate bar shape (the ‘barlense’ - lens¬ 
like components embedded in bars, Laurikainen et al.|2011 1 . 
In any case, we are not dealing with an extreme object. It 
is reasonable to expect that the conclusions above would 
be strengthened when considering more prominent bulges, 
while whether or not weaker ones would leave identifiable 
imprints may vary from case to case. It is not unlikely that 
for low-mass objects the differences induced in the kinemat¬ 
ics may fall within the grey area set by our ignorance on the 
B/P strength. 


6 FACE-ON SIGNATURES 


In this section we investigate the properties of the kinematic 
maps when the galaxy is viewed face-on. In this case the bar 
component is fully visible, while the B/P shape disappears 
in the projection. A recent comparisons between simulated 
and observed disk galaxies ( [Athanassoula et al.| |2014[ 
[Laurikainen et ah 20141 suggests that at low inclination 
B/P bulges turn into the aforementioned barlenses, opening 
a new route for face-on B/P diagnostics. These are crucial 
in order to establish the origin of B/P structures, the 
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Figure 17. End-on views of the disk component from simulations GTR106, GTRlll, GCS008 and CBL005. The first three were taken 
at snapshots such that their B/P strength was similar to that from GBL005 (disk+bulge) at t = 6 Gyr (see also Fig. |18| and |19[ |. Note 
that the FOV of GTR106 is larger than that of the other three runs. 



Figure 18. Major-axis kinematics for simulations GTR106, 
GTRlll, GGS008 and CBL005 in an end-on projection. The sim¬ 
ulations are taken at different times, as stated in the legend; these 
were selected so to have a similar B/P strength to CBL005 at 
t = 6 Gyr. The panels show the mean behaviour of the five kine¬ 
matic quantities along a 700-pc-wide slit centred on the 2 = 0 
plane. The results are derived from the disk component, except 
for GBL005 where both disk and bulge are considered together. 

properties of the galaxy, and to model correctly the density 
distribution along the line of sight (see, e.g.|Fragkoudi et al.| 

The reference work for the kinematic signatures of face-on 
B/P bulges is DOS. They investigate a number of possible 
B/P diagnostics in the density distribution and kinematics 
along the 2 direction and conclude that the most powerful 
proxy is the presence of negative double minima in hi. This 
result, based on a set of simulations of isolated disks, was 


Figure 19. As in Fig. |18| but for a slit centred on 2 = 1.5. 

later confirmed observationally by |Mendez-Abreu et al.| 

1 2008 ] ). 

Our set of simulations differs from that of DOS in sev¬ 
eral aspects (halo modelling, presence of gas and star 
formation - although both aspects are discussed in their 
paper - initial-condition setup) and we thought it in¬ 
teresting to compare our results to theirs. We bring an 
additional contribution to the study by looking at the 
evolution of this diagnostic with time for each of the runs 
and by investigating the behaviour of the stellar component. 

We discuss separately the strong and moderate-B/P 
groups of simulations, as well as the runs with star for¬ 
mation. In general, the behaviour along the major axis is 
reasonably representative of the whole maps; we therefore 
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discuss line plots mainly and show the 2D results only for 
a few representative snapshots. Panels showing the face-on 
maps for all the runs are available online (Fig. 10, 11 and 
12 for, respectively: the disk component at t = 0.1 Gyr, 
the disk component at t = 6 Gyr, the stellar component 
at t = 6 Gyr). As the odd moments are zero everywhere 
for symmetry reasons, we will only concentrate on the even 
ones (a and /14); the SB maps are visible online, while here 
we rely on the overlaid isodensity contours to hint to the 
visual appearance of the bar-B/P structure. In both the line 
plots and the maps, the filled symbols mark the location 
where the vertical extent of the structure reaches its maxi- 
mu «1!1 This is computed separately for the disk and stellar 
components. 


6.1 Simulations with strong B/Ps 


The ID results for this group of runs are shown in Fig. |20| 
One common feature is that until t ~ 3 Gyr the /14 profile 
is flat and consistent with zero. At t = 4 Gyr all the simu¬ 
lations have a B/P bulge, and the behaviour changes ac¬ 
cordingly. Negative regions develop in the centre, with broad 
and progressively deeper minima. Note how GTR106 tends 
to have the least striking changes with respect to the pre- 
B/P results; indeed, this is the simulation with the smallest 
B/P strength at these times. After t = 5 Gyr the evolution 
changes according to the specific history of each of the simu¬ 
lations. Both GTRlOl and GTR106 undergo a second rapid 
phase of B/P growth around t = 5.5 —6.5 and t = 6 — 6.5 Gyr 
(Fig. |2f and, correspondingly, the /14 curves start develop 
strong and clear minima roughly at the position of the B/P 
maximum-height location^ The importance of these fea¬ 
tures with respect to the initial, shallow minima can be 
fully appreciated when looking at the 2D maps (top row of 


Fig. 21 1 . These should be contrasted to those from GCS006 


(bottom row of Fig. 211; indeed, the B/P from this simu¬ 
lation does not experience multiple phases of rapid growth 
and the corresponding curves in Fig. [Revolve smoothly un¬ 
til the end of the simulation. Finally, GTR102 undergoes a 
second phase of growth only at the very end of the sim¬ 
ulated timespan; as a consequence, strong minima appear 
only at t = 10 Gyr and are only partially visible because 
they extend slightly beyond the FOV. 


6.2 Simulations with moderate B/Ps 

The ID results for this group of runs are shown in Fig. |22| 
In GTRlll, GTR116 and GGS008 the B/P shape starts be¬ 
ing visible between 5 and 6 Gyr, while in CBL005 within 


^ The vertical extent of the B/P is evaluated as median(|z|) and 
its maximum carries an intrinsic uncertainty of 150 pc, due to 
binning along the x direction. We compared this measure to the 
maximum of the 80th percentile of the \z\ distribution, as well as 
to the mean of the \z\ values between the 60th and 90th percentile, 
reaching good agreement. We do not show these two additional 
measures as they are indistinguishable from the first on the scale 
of the plots. 

® The featureless behaviour of GTRlOl at t = 10 Gyr (red curve 
at the bottom) is due to a temporary displacement of the minima 
away from the major axis. 
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Figure 20. Face-on behaviour of the disk component from the 
simulations developing a strong B/P over time. The panels show 
the mean value of along a 700-pc-wide slit centred on the 
y = 0 plane. The results are plotted at ten different times: from 1 
to 10 Gyr, in steps of 1 Gyr. The curves have been displaced 
vertically to facilitate their reading; the horizontal lines mark 
the = 0 level at each time and the tickmarks represent a 0.1 
interval in h^. Starting at t = 6 Gyr, the extent of the B/P is 
marked by a filled circle. 


the first gigayear. The minima deepens and extend out¬ 
wards as the B/P bulge grows. Their relative importance 
reflects the sequence in B/P strength (e.g. at late times: 
GTR116, GTRlll, GCS008, GBL005). As an example, the 
first two columns of Fig. |23| show the full maps for GTR116 
and GCS008 at t = 8 Gyr. 

For GBL005 Fig. shows the results from disk and bulge 
component taken together (solid lines), as well as the re¬ 
sults from the disk component only (dotted fine). As noted 
already by D05, the bulge tends to mask the minima; the 
result can be better appreciated by looking at the full 2D 
maps, because the masking effect is stronger away from the 
B/P major axis (third and fourth columns of Fig. 231. 


6.3 Stellar component 

The ID results for the stellar component of the runs with 
star formation are shown in Fig. |24| The first notable 
difference with respect to the two other cases shown 
above is in the presence of extended regions of positive 
/14, especially at early times. This is likely to be a side 
effect of the low velocity dispersion for this component 
and indeed a tight spatial correspondence between a and 
/14 features is generally present. At later times, regions 
of positive migrate outside the B/P area and central 
negative minima start developing. The behaviour of stars 
does not mirror entirely that of the disk component for 
the same runs. The second-buckling minima in GTR106 
develop earlier, are deeper and somehow better defined in 
the stellar component. Similar considerations hold for the 
other runs too, except for GGS008 where the h 4 minima are 
often well defined but not necessarily negative. They are 
still associated to the presence of a B/P bulge, though, as 
the two develop together starting at t = 6 Gyr. Overall, the 
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Figure 21. Full face-on kinematics of the disk component from simulations GTRlOl and GCS006 at t = 5,8 Gyr. The red, filled circle 
marks the B/P edges (see text). 



Figure 23. Full face-on kinematics of the disk component from simulations GTR116, GGS008 and GBL005 at t = 8 Gyr. For the latter 
also the version with disk and bulge together is shown. Meaning of the symbol as for Fig. |21| 
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Figure 22. Face-on behaviour of the disk component from the 
simulations developing a moderate B/P over time. Details as for 
Fig.[^ 


importance of the negative /14 minima follows once more 
the sequence in B/P strength: GTR106, GTR116, GTRlll 
and GCS008. 

Fig. 1^ shows the hi evolution of the two stellar popula¬ 
tions dehned in Sec. |4.4[ namely those born at t < 0.15 and 
t > 2 Gyr; the relative behaviour is similar for all the runs 
and features a progressive growth of positive regions at 
the centre and at the edges of the FOV. The most striking 
differences between the face-on kinematics of “old” and 
“young” stars are found for GTR116 at late times and we 
show their full 2D maps, together with those of the disk 
and full stellar component, in Fig. |26| To start with, we can 
appreciate the more prominent /14 minima characterising 
the stellar component with respect to the disk one - as 
discussed above. In general, the importance of hi minima 
is not affected by this or different age splits and indeed the 
feature is present in all stars, whether they form before, 
during or after the B/P. The central hi peak characteristic 
of the younger population is in spatial correspondence with 
a cr dip. 
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Figure 24. Face-on behaviour of the stellar component from all 
the simulations with star formation. Details as for Fig. |20| 
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Figure 25. Face-on behaviour of two separate stellar popu¬ 
lations; the one labelled as “old” comprises stars born at t < 
0.15 Gyr, while the “young” group concerns stars born at t > 
2 Gyr. The curves are shown starting at t = 3 Gyr. Other details 
as for Fig. 


6.4 Comments 

In summary, the main results of DOS for the disk component 
are confirmed, namely: (i) the presence of a B/P is accom¬ 
panied by negative /14 minima in the face-on kinematics 
and their importance is a function of the strength of the 
structure (ii) the B/P kinematic features are masked by the 
presence of a classical bulge. To this we add the kinematic 
consequences of a second buckling, i.e. the development 
of deeper, localised, external hi minima on top of central 
shallower ones related to the first phase of growth. We 
also inspect the behaviour of the stellar component, whose 
kinematics present some remarkable differences with respect 
to the disk. These are, most notably: (i) extended regions 
of positive /14 (ii) tendency to have deeper hi minima 
at a given time (iii) slight asynchrony in the evolution 
with respect to the disk component. We also note that in 



Figure 26. Face-on view of simulation GTR116 at t = 10 Gyr. 
The first row shows the disk component, while the bottom three 
different stellar samples: all stars, stars born at t < 0.15 Gyr and 
stars born at t > 2 Gyr, respectively. 


some cases stars may present B/P-related, non-negative hi 
minima. 

As for a, we showed the maps for completeness and 
performed a major-axis analysis similar to that presented 
for hi; in agreement with DOS, we did not find any signature 
that could point to the presence/extent of a B/P structure. 


7 EFFECT OF DIFFERENT SETUPS 

Simulations GTR102, GCS006 and GCS008 were added 
to the study in order to assess the effect of the choice in 
structural parameters - e.g. halo shape and density profile 
- on the results. GTR102 has a mildly non-axisymmetric, 
cored halo while GGS006 and GCS008 have a spherical 
NFW profile. They were picked, among their respective 
series, as they provide B/P histories similar to those from 
the GTR runs under consideration. Indeed, we do not want 
to assess the differences in the kinematics of objects with 
radically different evolutionary histories; we want to see 
if, for similar objects, differences in the underlying halo 
structure manifest themselves somehow. 

In addition to the halo differences, GGS006 and GCS008 
have a disk scalelength /i = 4 kpc, while the GTR runs have 
/i = 3 kpc. The GGSs also present a slower disk rotation 
than the GTR sample at the initial conditions. 


Notwithstanding the differences in the setup, the 
kinematic features turn out to be compatible when com¬ 
paring similar B/P shapes and strengths. Fig. 27 shows. 
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Figure 27. Disk component from simulations GTRlOl and GCS006 in a side-on and end-on view (top-two and bottom-two rows, 
respectively). Run GTRlOl is taken at t = 5 Gyr, while GCS006 at f = 8 Gyr; the choice is such that the two B/P structures present a 
similar morphology (boxy, in this case). 


as an example, the side-on and end-on views of GTRlOl 
and GCS006 at times where both simulations present a 
markedly boxy bulge. These maps can be compared to 
those presented for GTR102 in Fig. and |11[ as the run 
presents itself a boxy structure there. Note that although 
the B/P strength for GTRlOl is lower at the selected 
time, the morphology of the structures is akin among the 
three runs and so is the resulting kinematic behaviour. 
Similar considerations hold for the group of simulations 
with moderate B/P structures. Fig. 28 shows the case 
of GTRlll and GCS008 at a moment where their B/P 
bulges have compatible strengths. We show, in addition, the 
behaviour of the star particles alongside that of the disk. 
Again, the maps are very similar for both components. 

As for the face-on case, from what was shown and discussed 
in Sec. 1^ we can similarly conclude that, once the differences 
in the evolutionary history and strength of the B/P are 
taken into account, the kinematic features are compatible 
among the families of runs. 


Overall, appreciable differences in the behaviour of the 
kinematic moments - for similar structural properties - are 
only visible at early times. These mainly concern the impor¬ 
tance of the hs-V correlation, as already mentioned at the 
end of Sec. |4.1[ and are a result of the specific equilibrium 
reached in the initial conditions - whether it is the halo, the 
rotation properties of the disk or its density proHle. Flow- 
ever, these features soon leave room to the stronger ones 
induced by a B/P, even though, as discussed, some relics 
of the disk behaviour may survive in the global kinematic 
behaviour at small position angles and for weak B/Ps. Over¬ 
all, after their formation bars and B/Ps become the main 
drivers of the evolution of the kinematic maps and we do not 
see any evident trace in them that could lead to the initial 
structural properties of the disk/halo. 


8 THE EFFECT OF INCLINATION 


So far we have focussed on the line-of-sight kinematics for 
perfectly edge-on or face-on systems; we will now discuss 
how the B/P signatures evolve when progressively tilting 
the disk. This is particularly important for comparison 
to real data, for example samples of near-face-on galaxies 
such as DiskMass i Bershady et al.||2010 |. We will show the 
results from runs GTRlOl and GTR116 as representatives 
of the two groups of simulations under study and consider 
the disk component only, as the stars behave similarly. We 
chose, as starting points, the end-on and face-on case; we 
discuss the two scenarios in turns. 


The effect of inclination on major-axis, B/P diagnos¬ 
tics was studied by BA05 for a strong-B/P case in an edge- 
on projection; these were found to hold overall for i > 80 
{i = 90 corresponding to the edge-on case) and changes were 
explained in terms of the varying “screening” action of the 
disk at different projections. We recover their main hndings 
for the end-on case, namely: 

• the central density peak gains relative importance due 
to the progressive decrease of projected density in the outer 
parts; 

• the central velocity gradient becomes steeper, due to a 
lesser contribution of the (slow) disk to the LOSVD in the 
B/P region; 

• the central cr peak becomes broader, while the wings 
lose importance as the profile becomes dominated by B/P 
orbits; 

• the hs-V correlation turns into an anti-correlation as 
the bulk of the LOSVD moves to higher velocities; 

• the /i 4 minima become deeper. 

We find that these changes are induced very rapidly when 
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Figure 28. Simulations GTRlll and GCS008 in a side-on and end-on view (top-four and bottom-four rows, respectively). The behaviour 
of the disk and stellar components are shown separately. Run GTRlll is taken at t = 6 Gyr, while GGS008 at t = 7 Gyr; the choice is 
such that the two B/P structures have similar strength (based on the behaviour of the disk component). 


moving away from the end-on projection down to i ~ 75 
and then evolve rather slowly. In addition, we register 
variations in the importance of the central u peak (in 
opposite directions for GTRlOl and GTR116) and a 
considerable change in the hi shape. Especially some of the 
a results are better appreciated from a line plot; however, 
the quantitative assessment of the differences induced 
by variations in i is not what we are after here, and we 
prefer to show the behaviour of the full maps to grasp a 
few interesting 2D features. Though we have studied the 
variation of the maps in inclination steps of 5 degrees, we 
show only few rather inclined examples and comment on 
the results from all the other cases. 

Fig. |29| shows the behaviour of the runs at i = 90, 80, 70, 60 
and t = 6 GyiQ The major-axis changes outlined above 

® The case i = 90, i.e. the perfect edge-on views, were previously 


can be appreciated qualitatively even just from looking 
at the maps. A decrease in the central a value is visible 
in GTRlOl, while the opposite behaviour in GTR116 
escapes the eyes. The evolution in the hi shape is evident 
for GTRlOl and less striking for GTR116. In the former 
case, the edge-on morphology characteristic of strong B/Ps 
disappears for i < 85 as the minima collapse around the 
plane (although an X-shape survives until i « 70); in 
the latter, the minima simply shrink in size and become 
deeper. Interestingly, a rather considerable region of h^-V 
correlation survives away from the plane down to i ~ 70,80 
(for the strong and moderate-B/P case, respectively). The 

shown in Sec.[^ but we replicate them here to ease the comparison 
to the inclined cases. For the same reason we add the i = 0 case 
to the upcoming Fig. |30| even though this result was presented 
already in Sec.|^ 
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fourth and eighth rows of the panel show, as a sanity check, 
the behaviour of the respective disks at t — 0.05 Gyr, i.e. in 
the absence of any non-axisymmetric strncture. 

If we continue along the inclination sequence past i = 60, 
the maps evolve smoothly into those representing the 
face-on kinematics (see Fig. 21 23 or 30 as an example - 


note that the snapshots are different as well as the B/P 
position angle). Indeed, the trends in V, /13 and /14 slowly 
reverse: the central velocity gradient flattens, the strong 
hs-V anti-correlation regions progressively shrink and 
eventually disappear, while the central hi minima become 
shallower and shallower. 


D05 stndy the effect of inclination on the B/P diag¬ 
nostic for close-to-face-on disks. They conclude that the 
characteristic double minima in /14 can identify the presence 
of a B/P structure down to i « 30 (face-on being i = 0), 

shows the a and /14 maps^ 
for the two runs considered here at i = 0,10,20,30,40; 
the last row shows the result for the axisymmetric disk 
at i = 40 for reference. The a maps do not show ap¬ 
preciable evolution; this is confirmed by the line plots of 
the major-axis behaviour, except for a steady growth at 
the edges of the FOV which becomes even more evident 
at larger i. Regarding hi, in GTRlOl both the central 
(shallow) and external (deep) minima lose importance 
progressively. The former eventually disappear (around 
i ~ 55), while the latter evolve into those observed in the 
edge-on, side-on view (see e.g. Fig. albeit for a different 
snapshot). As for GTR116, the hi minima never really 
vanish either, although they similarly become shallower at 
f > 55 and evolve in shape towards those observed in the 
edge-on, side-on view (Fig.|^ again for a different snapshot). 

In Fig. and we showed also the kinematic 
behaviour of the disks at t = 0.05 Gyr to underline 
the contribution of the B/P structure at later times. At 
intermediate inclinations the disk is responsible for a bland 
h^-V anti-correlation at the most, while the hi maps are 
featureless. Bars and B/P bulges, conversely, are associated 
to strong /13 and hi features together with a richer a 
morphology. 


if the latter is strong. Fig. 


30 


However, as i decreases the characteristic B/P morphol¬ 
ogy is less and less visible and only the bar component can 


be clearly identified (see, however, Athanassoula & Beaton 
2006[ [Erwin fc Debattista|2013 1. 


This brings the question up again as to whether the presence 


of the B/P bulge can be guessed by the kinematics. Fig. 31 


shows the case of an end-on bar viewed at different inclina¬ 
tions near edge-on; we selected GTRlOl at t = 2.5 Gyr as 
it provides one of the B/P-less cases with the strongest bar, 
but equivalent considerations hold when we inspect bar-only 
snapshots from GTR116. Overall, we find extended similar¬ 
ities in the line-of-sight kinematics for the two cases - with 


We do not to show the V and hg maps in Figs, [soj and [ 3 ^ - 
namely those starting from a face-on projection, where the two 
quantities are null by definition - to avoid overcrowding. The 
evolution of these maps from the face-on case to more inclined 
projection is gradual and progressively leads to the features dis¬ 
cussed in Sec. [ 4 ] for the edge-on case. 


and without a B/P. Down to inclinations i « 80, the ap¬ 
pearance of the surface-brightness map can still discriminate 
between the two even though the other kinematic quantities 
behave likewise. We stress once more that if the spheroidal 
appearance of the isodensity contours originated from a clas¬ 
sical bulge, then the kinematics would be rather different; 
the presence of a B/P is therefore easy to establish in this 
case. At lower i, the bar-only and B/P maps become yet 
more similar and the morphological criterion less effective 
in distinguishing the two. The shape of the hi minima is 
such that these avoid the centre in the B/P case and in¬ 
stead fill entirely the contours in the bar-only case - some¬ 
thing which becomes even more evident at i values lower 
than what shown - but overall the signatures resemble each 
other. It therefore may not be straightforward to identify a 
B/P on top of an end-on bar at inclinations i < 80. 
However, in the side-on case the situation is different. We do 
not show and comment explicitly the side-on case at nearly 
edge-on inclinations, because the characteristic B/P shape 
is still somewhat visible there - especially in GTRlOl. Any¬ 
ways, besides the surface brightness maps, the importance 
of the hi minima is still rather different in the B/P and 
B/P-less cases at those inclinations. Instead, we move on to 
nearly face-on projections. Fig. already showed how the 
characteristic hi minima associated to a face-on B/P resist 
the effect of disk tilting; Fig.j^now presents the same maps 
for the bar-only case. There is no doubt that a bar alone has 
a radically different line-of-sight kinematics at these projec¬ 
tions. This holds at all i values down to the edge-on, side-on 
case; we can therefore conclude that when viewed side-on 
bars and B/Ps have different signatures regardless of i. 

It is outside of the scope of this paper to carry out a full 
inclination study; we remark, however, that the above con¬ 
clusion is consistent with what found by D05, namely that 
a position angle < 45 is more favourable for finding B/P 
bulges in inclined systems. 


9 SUMMARY AND CONCLUSIONS 


In this paper we showed the 2D, projected line-of-sight 
kinematics of bars and boxy/peanut bulges from a set 
of recent dynamical simulations featuring star formation. 
An average of a thousand Voronoi-binned kinematic maps 
were produced for each of the runs, in order to capture the 
time evolution for at least eight different viewing angles 
(face-on, edge-on with the bar/peanut at seven different 
position angles). The information contained in this dataset 
exceeds what can be discussed in a single paper and here 
we presented an initial selection of the results, building 
on previous works studying the kinematic signatures of 
boxy/peanut bulges. 


Bureau & Athanassoula (20051 identified a set of character¬ 


istic imprints induced by simulated boxy/peanut structures 
on the ID, major-axis kinematics of edge-on galaxies - 
primarily the presence of a correlation between the first 
and third moment, /13 and V. We recover their results 
in both the disk and stellar component, as well as their 
dependence on the strength of the boxy/peanut and its 
position angle. To this we add the identification of strictly 
peanut-related signatures in rough spatial correspondence 
with the projected edges of the structure, away from the 
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Figure 29. Effect of disk inclination for simulations GTRlOl and GTR116 (top-five and bottom-five rows, respectively). The starting 
point of the study is the end-on view at t = 6 Gyr. The fifth and tenth rows show the behaviour of B/P-less disk at t = 0.05 Gyr and 
for an inclination of 60 degrees. 


major axis. These are in the form of elongated wings of 
large hs values and X-shaped regions of deep minima. 
Both become more evident at large boxy/peanut strengths; 
the former is mainly characteristic of intermediate position 
angles (30 — 60 degrees, roughly), while the latter are 
visible down to the end-on view. The extent to which these 
features can be used to constrain the position angle of the 
bar/peanut system is limited and varies with the strength 


of the latter. For prominent structures, the amplitude and 
morphology of ha and hi features undergo a considerable 
evolution when rotating the system from a side-on to an 
end-on view and this, in turns, should allow the bracketing 
of the position angle in few rough categories. For moderate 
boxy/peanuts, however, the evolution with position angle is 
slow and does not feature outstanding signatures, making 
it more difficult to draw conclusions on the orientation of 
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Figure 30. Effect of disk inclination for simulations GTRlOl and GTR116 (first two and last two columns, respectively). The starting 
point of the study is the face-on view at t = 8 Gyr. The last row shows the behaviour of B/P-less disk at t = 0.05 Gyr and for an 
inclination of 40 degrees. 


the structure. 

In the end-on view, the projected morphology of a 
boxy/peanut resembles that of a classical bulge and, 
notwithstanding the expected differences in the light 
profiles, discerning between the two structures may not 
be straightforward. Here we investigate the case where 
both are present and how the projected kinematics of 
the boxy/peanut is altered by the presence of a classical, 
spherically-symmetric component. We find that the char¬ 
acteristic signatures are considerably weakened and would 
hint to a mnch more moderate structure than what really 
present; that the observed morphology cannot be generated 
solely by a peanut is however evident when comparing 
the kinematic maps to cases with a similar projected 
appearance and no classical bulge. 

In the face-on case, |Debattista et al.|p005[ ) had identified a 
peanut diagnostic in the presence of symmetric, negative /14 


minima. We confirm this result and witness the development 
of even deeper, narrower minima spatially confined around 
the peanut edges; these occur in the simulations with strong 
boxy/peanut structures if and when they undergo a second 
buckling phase. 

The analysis of the disk component made the backbone 
of the paper and results for the stars formed during the 
simulations were discussed as an aside. The reason for 
this is in that the two components present globally similar 
behaviours concerning their line-of-sight kinematics. Stars, 
however, tend to carry more extreme amplitudes for the 
higher moments even if the overall morphology of the maps 
is never too dissimilar from that of the disk; this is true 
in both the edge-on and face-on views. In addition, stars 
born at radically different times during the simulation 
(e.g. before and during/after bar and peanut formation, as 
was shown here) may present, on top of a different spatial 
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Figure 31. Effect of disk inclination for simulation GTRlOl at t = 2.5 Gyr, i.e. when only a bar is present. This is viewed end-on in 
the top row. 



Figure 32. Effect of disk inclination for simulation GTRlOl at 
t = 2.5 Gyr, i.e. when only a bar is present. This is viewed face-on 
in the top row. 


distribution, also distinguished kinematic behaviours. This 
aspect will be investigated further elsewhere. 

Perfectly edge-on and face-on views are just the two 
extremes of all orientations a disk galaxy can possibly 
be viewed at. It is therefore desirable to understand how 
the signatures discussed above evolve as the disk is tilted 
away from these two projections. As already remarked by 
Debattista et al. ( 2005[), in inclined systems boxy/peanut 


bulges are more difficult to identify at large rather than at 
small position angles. In agreement with this, we find that 
in the end-on view the signatures become soon very similar 
to those of a peanut-less bar when i < 80. Conversely, 
however, when a side-on boxy/peanut is inclined away from 
a perfectly face-on view, the characteristic /i 4 minima sur¬ 
vive at each resulting projection - albeit with progressively 
lower strengths. 

The bulk of the simulations we discussed comes from the 
GTR series of Athanassoula et al. (20131 with a cored, 
spherical dark-matter halo. These simulations differ from 
one another just in the initial mass fraction of the disk which 
is in the form of potentially star-forming, gas particles. In 
order to assess the scatter induced in the results by different 
system setups, we added a few different simulations to the 
study: another GTR simulation characterised by a mildly 
non-axisymmetric halo; two simulations with a cuspy halo 
and different disk scalelengths; a case with a classical bulge, 
mainly serving the purpose of assessing its masking effect. 
Although we do find differences in the resulting kinematics 
that can be traced back to the different initial rotation 
curves, these are only really appreciable during the pre¬ 
bar/peanut stages; once the latter form, they become the 
main drivers of the behaviour of the kinematic moments. 
Therefore, even though halo properties are known to play 
a major role in determining the evolution of bars and 
peanuts, similar structures will give rise to compatible maps 
- regardless of the features of their respective halo and disk. 
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The only exception to this is in the amplitude (but not 
the morphology) reached by o, which retains memory of 
the different initial rotation properties of the disk. Also in 
the relatively bland, side-on maps for the highest moments 
we register some small case-to-case variations that may 
be traced back to different initial setups. These concern 
mainly /13 and are, however, without major consequences 
besides complicating further the position-angle diagnostics 
for moderate boxy/peanuts. 

In summary, our results can be used to better identify 
bars and boxy/peanut structures from 2D kinematic data 
as those provided by latest generation IFUs. Indeed, these 
structures leave distinct photometrical imprints only in a 
limited range of viewing angles and it may be difficult to 
clearly pinpoint their presence - let alone their position an¬ 
gle and other properties - from the light distribution of the 
object. The strength of dynamical simulations in accurately 
following the formation of boxy/peanut bulges is traded, of 
course, with the lack of a realistic formation history of the 
simulated object. We therefore stress that our results de¬ 
scribe an idealised and unperturbed setting; an assessment 
of how they would change if all the processes relevant to 
galaxy formation were taken into account is currently be¬ 
yond our reach. 
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